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0\ . We assess the feasibility of probing dark energy with strong gravitational lensing systems. 

'/^ ' The capability of the method, which depends on the accuracy with which the lensing systems 

I are modeled, is quantitatively investigated using the Fisher-matrix formalism. We show 

that this method might place useful constraints on the density parameter and the redshift 
evolution of the dark energy by combining it with a constraint from supernova measurements. 
For this purpose, the lens potential needs to be precisely reconstructed. We determine 
fH \ the required quality of data. We also briefly discuss the optimal strategy to constrain the 

^ , cosmological parameters using gravitational lensing systems. 

o , 

. §1. Introduction 

Recent observations of distant Tvpe la supernovae have provided strong evidence 
for the acceleration of the universe. cP'B) Recent measurements of cosmic microwaYe 
^ ' background anisotropics favor a spatially flat universe with cold dark matter, B^'cP 

. while measurements of the large-scale structure in the distribution of galaxies favor 

a low-density universe. These observations can be explained by the hypothesis 
that our universe is dominated by dark energy in addition to cold dark matter. 
Motivated by these findings, many authors have recently investigated the origin of 
dark energy. i^'B'S* An attractive feature of these theoretical models is that the 
coincidence problem, the near coincidencfi^ of the density of matter and the dark 
mergy component may be explained, EP'@)'t3) though the plausibility is under debate. 



An evolving scalar field has beeri investigated as a model of a decaying cosmo- 
logical constant (dark energ energy can be characterized by the 
effective equation of state w = pd/J^d-, where po and Pd are the energy density and 
the pressure of the dark energy, respectively. Different theoretical scenarios of dark 
energy may be distinguished by the measurement of the effective cosmic equation 
of state. For example, the quintessence model generally predicts that w can be a 
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function of redshift and satisfies w > —1.0) If darkjjnergy originated from the cosmo- 
logical constant or a false vacuum energy, w = — 1. lP Parker and Raval have explained 
dark energy in terms of a quantum vacuum effect produced by the curved spacetime 
of general relativistic cosmology. In their model, u; < — 1 is allowed. Models which 
predict ■(i2< —1 have also been proposed within the framework of the evolving scalar 
field. Therefore, the measurement of w may provide information regarding the 

origin of dark energy. For this reason various observatiorial stiat^ies have been pro- 
posed to probe characteristics of dark energy. lla^'Ej)'LH)'ElP'Bj)'E3)'L3^'c3''Ej)'E3^ One 
of the most promising approaches for determining w employs type la supernovae, 
whik other approaches using CMB anisotropies and number-counts should be useful 
too.Ba^'CJ^'Ey) However, it has been pointed out that the method using supernovae 
is fundamentally limited because the luminosity distance depends on w through a 
multiple integration, which smears out information concerning w and its time vari- 
ation. c3 

With the situation as described above, it is very useful to have an independent 
method to determine w. One possible such method is to use gravitational lensing. In 
fact, it is known that lensing statistics depend strongly on the value of the cosmolog- 
ical constant, as well as the characteristics of dark energyJ^iS^'a^'l^ Another 
possible method to probe dark energy has been proposed. E^^'u^'S'eS This method 
relies on precise measurements of 'clean' gravitational lensing systems. Here 'clean' 
means suitable for modeling a gravitational lensing system, as is an Einstein ring 
or Einstein cross. Motivated by this recent proposition, we assess the feasibility of 
the method to probe dark energy that employs the Fisher-matrix formalism, which 
is useful to demonstrate how accurately one can estimate model parameters from a 
given data set. (For a review see, e.g., Ref.^).) The primary goal of the present 
paper is to determine the quality necessary for a data set from gravitational lensing 
systems in order to probe dark energy. 

A gravitational lensing system can be used to measure the ratio of (angular 
diameter) distances, while the supernovae can be used to determine the luminosity 
distance itself. One might ask whether the method of the gravitational lens provides 
new independent information. The second goal of the present paper is to answer 
this question. We show that the gravitational lens method can be a useful tool that 
complements the supernova method as a probe of dark energy. 

This paper is organized as follows. In §2, we briefly review basic formulas. 
Then, we investigate the capability of the method employing strong gravitational 
lensing systems, focusing on errors in determining cosmological parameters, in §3. 
Uncertainties involved in modeling a lens galaxy are discussed in §4. In §5, we 
discuss the redshift distribution of a lensing system for the purpose of optimizing 
the method. Section 6 is devoted to a summary and conclusions. Throughout this 
paper we use units in which the velocity of light, c, equals 1. 

§2. Basic formulas 



In this section we give a brief review of basic formulas for the topics considered 
in the present paper. We restrict ourselves to a spatially flat FRW universe, in which 



3 



case the Friedman equation is written H{z)^ = i/g [■!^o(l + z)^ + (1 ~ ^o)f{z)]j where 
Hq = 100/ikm/s/Mpc is the Hubble constant, ]7o is the cosmic density parameter of 
the matter component, and f(z) describes the redshift-evolution of the dark energy. 
When the cosmic equation of state is parameterized as w{z) = wo(l + z)'^ , we have 



f{z) = (l + zfexp 



3wo 



'{l + zY-l 



(2-1) 



which reduces to f{z) = {1 + z)'^^^~^'^°^ in the hmit — > 0. In a spatially flat universe, 
the angular diameter distance between zi and Z2 is 

1 dz' 

= Ho{l + Z2) A, [noil + z'r + {l-Qo)f{z')V/^' ^^'^^ 

where we have assumed zi < Z2- 

In order to estimate the accuracy to which we can constrain Qq, wq, and 
with measurements using gravitational lensing systems, we employ the Fisher- matrix 
approach. With the Fisher-matrix analysis, one can estimate the best statistical 
errors on parameters from a given data set. For this reason, this approach is widely 
used to estimate how accurately cosmological parameters are determined from the 

arge scale structure of galaxies or the cosmic microwave background anisotropics. 

ii^'E^'C^P Tegmark et al. applied the Fisher-matrix analysis to supernova data sets. 

!S Their analysis is useful for the present study. In general, the Fisher-matrix is 
defined by 

where L is the probability distribution function of a data set, given model parameters 

ei. 

To compute the Fisher matrix, we need to make an assumption concerning the 
probability function L, and more specifically concerning the statistics of errors of 
measurements. The method employing a strong gravitational lensing system relies on 
the assumption that the ratio p£the angular diameter distances, Dls/Ds, is precisely 
determined for each system, E3)'E3) where Dls is the angular diameter distances 
between the lens and a source object, and Ds is the distance between the source 
and the observer. A simple example in which Dis/Ds is determined with a lensing 
system is considered in §4. Here, we assume that the lens model can be determined 
by measurements of images of a lensed object and the velocity dispersion of the lens 
galaxy, as well as the redshifts of the source and the lens objects. In a realistic 
situation, observational data is contaminated with errors in modeling the lensing 
potential and measurement of the velocity dispersion, as we discuss in §4. Hence, an 
observed value of Dls/Ds contains error. Throughout this paper, we also assume 
that this error can be described by a Gaussian random variable, for simplicity. For 
definiteness we assume that a data set consists of N gravitational lensing systems 
and that the observed ratio of the angular diameter distances of the n-th lens system 
can be expressed as Rn = Disn/ Dsn + ^n, where Disn and Dsn are the angular 
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diameter distances Dls and -^^5 of the n-th lensing system, and the error Z\„ behaves 
as a Gaussian variable with zero mean, {An) = 0. In this case the Fisher matrix 
(|2-3|) can be written 



N 



1 



n=l\ "/ 



1 



N 



(2-4) 



n=l 



where Wni = 91n (^Disn/Dsn^/dOi, En describes the variance of the Gaussian vari- 
able with (A^) = {DiSn/ Dsn)^£n^ the last approximate equality holds in the 
case that £«(= e) <C 1 for all n. In contrast to the case of supernovae,o^ the fac- 
tor 2 appears on the right-hand side of Eq. (2-4). The corresponding term can be 
ignored in the case e <^ 1. This factor is traced back to the assumption that the 
error Z\„ is proportional to Disn/ Dsn- By using the Bayse theorem, the probability 
distribution in the parameter space can be written 



P{9i) (X exp 



Y.{0i-et)F,,{9,-e^ 



(2-5) 
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Fig. 1. Confidence regions from various data sets for gravitational lensing systems. Here we have 
chosen the target parameters as — 0.3 and wq — —1, with the number and the individual 
uncertainties as follows: (a) N/e^ = 10^ (b) N/e^ = 3 x 10^ (c) N/e'^ = 10*; (d) N/e = 3 x 10*, 
where we have set e„ equal for all systems, £„ — e. The curves are the contours of the 68%, 
95%, and 99.7% confidence regions. We chose the redshift distribution of lensing systems zg and 
2^ through a homogeneous random process in the ranges 1 < zg < 3 and 0.3 < zi < zg — 0.5, 
respectively. 
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where we have assumed that errors in the target model parameters Oj^ are small. 
Thus the Fisher matrix gives the uncertainties in the parameter spaces, which are 
described by a Gaussian distribution function around Oj"^. 

§3. Estimating errors in determining cosmological model parameters 

We now demonstrate how accurately cosmological model parameters can be de- 
termined by using the Fisher-matrix formalism. We first consider the model in which 
the cosmic equation of state is constant {w{z) = wq), which is equivalent to the as- 
sumption f{z) = {l + z)^^^'^'^°\ Considering the Fisher-matrix to be the 2x2 matrix 
corresponding to and wq, we show how accurately these two parameters can be 
determined. Figure 1 displays the confidence regions in the f2o-wo plane described by 
the probability function (2-5), where we used the cosmological parameters i7o = 0.3 
and Wq = —1 as the target model parameters. We considered various data sets for 
the number and the variance of the errors in the gravitational lensing systems. The 
quality of the data is characterized by N/e'^, because the Fisher-matrix scales in 
proportion to N/e^, whose value is shown on each panel of the figure. For the red- 
shift distribution of the lensing systems, we chose zs and zl through a homogeneous 
random process in the ranges 1 < zs < 3 and 0.3 < zl < zs — 0-5. The three curves 
in each panel indicate the 68%, 95%, and 99.7% confidence regions. Figure 1 shows 
that the method of the gravitational lens is sensitive to J?o, but rather insensitive 
to Wq in the case that N/e"^ is small. It also shows that a useful constraint on f2o 
and Wq can be obtained if many lensing systems are measured precisely and N/e'^ is 
sufficiently large. 

One might think that the gravitational lens method can be no better than the 
supernova method because the latter supernova measures the luminosity distance 




Fig. 2. Confidence regions from a data set from supernovae with individual uncertainties of 0.15 
magnitude. Here, 100 supernovae randomly distributed in the range 0.5 < 2; < 1.5 are consid- 
ered. The curves are the contours of 68%, 95%, and 99.7% confidence regions, and the same 
target model as Fig. 1 is used. 
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Fig. 3. Same as Fig. 1, but used combining the constraint from the data set of 100 supernovae 
displayed in Fig. 2. 



itself, while the former measures only the ratio of the (angular diameter) distances. 
However, in principle, these two methods give independent information and can be 
complementary. Figure 2 demonstrates how accurately we can put a constraint on 
model parameters using a data set from 100 supernovae randomly chosen in the range 
0.5 < z < 1.5 with individual statistical uncertainties of 0.15 mag.C3^ By comparing 
Figs. 1 and 2, it is clear that the two methods provide independent constraints. Thus, 
using both methods, we can eliminate ambiguities in determining model parameters 
that arise when using only one. Figure 3 displays the confidence regions in the Qq- 
wq plane, which is the same as Figure 1, but with the addition of the constraints 
obtained from the data set of 100 supernovae presented in Fig. 2. It is clear from 
Fig. 3 that the ability to constrain and wq is significantly improved. 

We next consider the case that the density parameter Qq is fixed from other 
experiments. We assume i7o = 0.3. By considering the Fisher matrix to be the 2x2 
matrix corresponding to the parameters wq and v, Fig. 4 displays the confidence 
regions on the wq-u plane, where we used wq = —1 and = as the target model 
parameters. This figure displays the results using both the gravitational lensing 
systems and the data set of 100 supernovae, for which we assumed the same range 
of distribution and uncertainty of sources as in Fig. 2. The constraint on v might 
not be strict when N/e'^ is of order 10^, but the useful constraint is obtained when 
N/e'^ ~ O(IO^). Hence precise measurements of many gravitational lensing systems 
might provide a useful probe, even for the parameter v. 
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Fig. 4. Confidence regions in the wo-i' plane, with density parameter fixed as ilo ~ 0.3. The target 
model parameters wq = ~1 and v — are used. The situation regarding the redshift distribution 
and N/e^ are the same as in Fig. 1. 



§4. Errors in modeling lens galaxies 

In this section we discuss e, which characterizes systematic errors in the ratio 
Dis/Ds- For this purpose, we start by demonstrating one of the simplest ways 
to determine the ratio, adopting a lens potential based on an isothermal ellipsoid 
model. H'E§' In this model, the lens equation gives an elliptical image of the Einstein 
ring with minor and major axes 

0± = OeVi^I, (4-1) 

with 

Us 

where is the one-dimensional velocity dispersion, and the ratio e of the minor axis 
to the major axis is related to the ellipticity e by e = ^(1 + e)/(l — e). Thus the 
ratio Dls/Ds can be determined by measurements of e, 9e and dy. 

We note that several authors haverecently pointed out possible errors involved 
in modeling the lensing potential. e3)'e3)'E§)'e£P With regard to this matter, we first 
consider a possible density profile that differs from the singular isothermal profile, 
though we restrict ourselves to a saherically symmetric density profile. We assume 
a generalized Navarro-Frenk-WhitetJ^'ElP or Zhao density profile,! 



52h 



Pir) = Vf^, (4-3) 
g{r/rs) 
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where we have defined 

g{x) =x"{l + x)^-'^, (4-4) 



where a and are parameters and Vg is a characteristic length scale (cf. ref. ^ 
In this case, the Einstein ring radius is determined by solving the lens equation 



2DlS 



eE-Fi9E,DL/rs,a)i7ra',-^ = 0, (4-5) 



where we have defined 



1 r^E 1 r°° Dr/r 
F{eE,DL/rs,a) = — dO'e'- dt , ' ^. (4-6) 



Oe Jo vr J-oo giyW^^T^Di/r, 



Here we have considered the case C = 2, so that the lens model possesses the same 
velocity dispersion as the singular isothermal sphere at large distances. In the limit 
a — > 2, we have F{6e, Di/rs,a) = 1, and the model reduces to the singular isother- 
mal model. This model indicates how the ratio Dis/F)s, which is determined by 
solving the lens equation (4-5), depends on the density profile parameterized by 
and a with the same observable quantities 9e and o"^. Figure 5 plots the normalized 
difi'erence between Dls/Ds and the corresponding value in the singular isothermal 
model. Explicitly, this figure displays contours of the quantity 

Dls/Ds -Dls/Ds\sis ^ 1 ^ 

Dls/Ds\sis F{eE,DL/rs,a) ' ^ ' ' 

in the r^-a plane, where Dis / Ds\sis = ^£;/(4vr(T^) is the value of this ratio in the 
singular isothermal model. In this figure we fixed Be = ^" and used the cosmological 
model with density parameters Hq = 0.3 and w{z) = —1. The four panels correspond 
to the cases in which the redshift of the lens is zl = 0.5, 1.0, 1.5 and 2.0. For 
example, the normalized difference is of order 5% and 10%, respectively, for a = 
1 and a = at = 0.1 /i~^kpc, which depends weakly on the redshift of the 
lens objects. The difference between the ratios becomes larger as increases or a 
decreases. We thus find that the ambiguity in the density profile can be an important 
factor in modeling the lensing potential. 

We next consider the errors involved in the observation of the velocity dispersion, 
which we have not yet discussed. Mortlock and Webster systematically investigated 
the errors involved in gravitational lensing statistics and pointed out that the greatest 
uncertainty comes from the dynamical normalization of lens galaxies. Ej* According 
to their results, the velocity dispersion in the isothermal sphere model (the line-of- 
sight velocity dispersion away from core radius) can be related with the observed line- 
of-sight velocity dispersion of stars cT||(i2/), assuming a circular aperture of radius 
Rf as 

a, = a\\{Rf)(A + B^\ (4-8) 



9 




Fig. 5. Contours of the normalized difference defined by Eq. (4.7) in the Ts-ot plane. This figure 
displays the normalized difference between the ratio Dig/Dg and the corresponding value for 
the singular isothermal model. Note that a = 2 represents the case of the singular isothermal 
model. We fixed 6e ~ 1" and adopted the cosmological model with Hq — 0.3 and w{z) = —1. In 
each panel, the redshift of the lens object is adopted as = 0.5, 1.0, 1.5 and 2.0, respectively. 
The curves are the contours of the 1%, 2%, 5% and 10% levels from left to right for all panels. 

where A and B are constants independent of the core radius Tc, axiii Rg is the 
effective radius of the lens galaxy. Kochanek used A = 1 and B = 2.E3) Then, the 
non-dimensional uncertainty of may be expressed as 

a^iiRf) ^ A + BrjRg' ^ ' 

where 5a\\{Rf) and 5rc represent uncertainties on a\\{Rf) and r^. Because the ratio 
of the distance Dis/Ds is proportional to cr^, the error that contributes to e is 
estimated as e ~ 2(5(T^ /cr^, where we have assumed ba^joy <^ 1. Concerning the 
uncertainty of the steller velocity dispersion, measurements of the velocity dispersion 
of lensing galaxies with acciijacY to within approximately 3% have been made with 
the Keck-II 10m telescope, E^^'l^ while the uncertainty on the core radius rc might 
be very important. For example, the second term on the right-hand side of ( [4 -91 ) is 
of order 0.1 for = 0.1 kpc and Rg = 2 kpc. Thus the investigation of this section 
demonstrates the importance of precisely modeling the lensing galaxy, taking the 
possible finite core radius and the density profile into account in a realistic analysis. 
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§5. Optimizing constraint on model parameters 

In this section we address the problem of determining the optimal distribution 
in redshift of gravitational lensing systems in order to best constrain the cosmo- 
logical parameters. This problem might be only of theoretical interest, but having 
such information could be useful in planning observation of gravitational lensing sys- 
tems. The same problem has been considered by Huterer and Turner with regard 
to supernova data.Ea' Our discussion here is based on the following observation: By 
a theorem concerning the Fisher matrix, it can be shown that the minimum error 
attainable on 6i is expressed by the diagonal part of the inverse Fisher-matrix if the 
other parameters are known. clP To address the above stated problem, we investigate 
what the pair of redshifts of the source and the lens objects maximize the diagonal 
part 

where the case of one lensing system is assumed for simplicity. It is useful to note 
how the normalized variation of the ratio Dis/Ds is related with the variances of 
the cosmological parameters and the Fisher-matrix: 

^g^^ = E W.Ae. = sFj^ilAno + sF^l^Awo + eF^'Au. (5-2) 

1/2 

Figure 6 displays the amplitude of the diagonal part given in eF^^' as a function 
of zl for zs = 5, 4, 3, 2, 1 and 0.5. The left and right panels depict the cases of the 
target model with i7o = 0.3, wq = —1 and = and with fig = 0.3, wq = —1 and 
v = —1/2, respectively. The top, middle, and bottom panels, respectively, display the 
components 9i = Qq, wq, and u. The amplitude shown in each panel confirms that 
l7o can be most easily determined, while wq and are relatively difficult to determine 
than l7o- From the top panels, it is clear that i7o is more precisely determined as zl 
becomes larger and that the precision with which l7o can be determined is almost the 
same for all zl satisfying zs > zi }^ 1. Contrastingly, the situation for wq and u is 
complicated. We see that as the redshifts of the source and the lens objects become 
larger (e.g., for zs > zl ^ 2), the precision with which wq and v can be determined 
becomes better. However, this is not always the case, because Fa becomes zero for 
certain combinations of zl and zs, which depend on the target model parameters. 
From the middle and lower panels, we can find another peak for zs ^ I at which 
Wq and u can be determined precisely. It is thus seen that low redshift lensing 
systems are also useful to determine the redshift evolution of dark energy. This is 



consistent with the result of Ref. 39), and it is also consistent with the result obtained 
from supernova analysis. £3-' This feature can be understood from the fact that the 
dark energy component increases relative to the matter component as the redshift 
decreases. 
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1/2 

Fig. 6. Amplitude of eFJ as a function of the redshift of the lens object 2^ with various redshifts of 
the source object zg- The top, middle and bottom panels, respectively, display the components 
Qo, Wo and v. The left and right panels correspond to the case of the target model with i?o = 0.3, 
Wo = —1 and v = Q and with Qq = 0.3, wo = —1 and v = —0.5, respectively. In each panel, the 
curves correspond to zg = 0.5, 1, 2, 3, 4 and 5, from top to bottom near zi, = 0. 



§6. Summary and conclusions 

In summary, we have assessed the feasibihty of probing dark energy using gravi- 
tational lensing systems. Wc have demonstrated how precisely the model parameters 
can be constrained using the Fisher-matrix formalism. Our results show that the 
method might place useful constraints on the density parameters and the cosmic 
equation of state if many lensing systems are measured precisely. Our results also 
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demonstrate the toughness of the method. It was found that N/e^ ^ 10^ is required 
to probe the cosmic equation of state wq even if combined with the supernova method 
employing data from 100 supernovae, and N/e'^ ^ 10^ is required to determine wq 
using only gravitational lensing systems. 

In order to realize the condition N/e'^ = 10^ - 10"^ we need e ^ 0.1 - 0.3 for 
N = 100. This required accuracy of measurement and modeling the lensing sys tem 
might be difficult to realize, due to theoretical and observational ambiguities. £5* As 
discussed in §4, uncertainties involving the density profile and the velocity dispersion 
are problematic for precise modeling of a lensing system. Nevertheless, the possibil- 
ity of precisely modeling lens potentials has been studied. For example, the ability of 
observed Einstein rings to precisely model lens potentials has been discussed by sev- 
eral authors. eH^'E^P The advantage of Einstein rings is the richness of the information 
they provide. The authors of those works have claimed that the distribution of the 
magnifications of Einstein rings makes it possible to eliminate lens model ambiguities 
and to model the mass distribution of a lens precisely. S Furthermore, they have 
pointed out that Einstein rings must be found frequently in multiply imaged quasars 
with sufficiently long observations, because it is likely that all quasars and AGN have 
host galaxies. Thus we might expect that studies of modeling lensing galaxies will 
provide more accurate results in the near future, employing many sample galaxies 
in SDSS, and that future progress might allow us to solve the problems involved in 
modeling the lensing potential. 

In §5, we considered optimizing the analysis. Though our consideration was 
restricted to the case of one gravitational lensing system, it is worth comment. We 
found that the gravitational lensing system method is most sensitive to l7o, and less 
sensitive to wq and i^. To determine Hq, it is most advantageous to have larger 
redshifts of the source and the lens objects, while this is not always true for wq and 
z/. This feature is due to the fact that the dark energy component decreases relative 
to the matter component as the redshift increases. Lensing systems with smaller 
redshifts can be useful to probe dark energy. 
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